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ABSTRACT 

The nature of the ultra-luminous X-ray sources (ULXs) in the nearby galaxies is a matter 
of debates. One of the popular hypothesis associates them with accretion at a sub-Eddington 
rate on to intermediate mass black holes. Another possibility is a stellar-mass black hole in 
a high-mass X-ray binary accreting at super-Eddington rates. In this paper we find a highly 
significant association between brightest X-ray sources in the Antennae galaxies and stellar 
clusters. On the other hand, we show that most of the X-ray sources are located outside of these 
clusters. We study clusters associated with the ULXs using the ESO Very Large Telescope 
spectra and the Hubble Space Telescope data together with the theoretical evolutionary tracks 
and determine their ages to be below 5 Myr. This implies that the ULX progenitor masses 
certainly exceed 40 and for some objects are closer to 100 solar masses. We also estimate 
the ages of clusters situated close to the less luminous X-ray sources (with luminosity in the 
range 3 x 10 38 < Lx % 10 39 erg s _1 ) and find that most of them are younger than 10 Myr, 
because they are surrounded by strong Ha emission. These findings are consistent with the 
idea that majority of ULXs are massive X-ray binaries that have been ejected in the process 
of formation of stellar clusters and at the same time rules out the proposal that most of the 
ULXs are intermediate mass black holes. The very small age of the clusters associated with 
ULXs favours the dynamical few-body encounters in the dense cores of forming star clusters 
as the main ejection mechanism. 

Key words: accretion, accretion discs - galaxies: individual (NGC 4038/NGC4039) - galax- 
ies: star clusters: general - X-rays: galaxies 



1 INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are non-nuclear X-ray 
sources with luminosities exceeding an Eddington limit for a lOM© 
black hole (see Feng & Soria 201 1 for a review). The most popular 
models for the ULXs involve either intermediate mass black holes 
(IMBH) with masses of 10 3 - 1O 5 M or stellar -mass black holes 
(StMBH) accreting at highly super-Eddington rates. Both scenar- 
ios require massive donors in a close orbit (similar to SS 433, see 
Fabrika 2004 ), and the later model also requires some collimation 
of radiation into our line of sight ( Fabrika & Mesc heryakov|200T| 
King et al. 2001 ; Poutane n et al.|2007| ) to explain the objects up to 
possibly 10 4i erg s _i . 

The X-ray luminosity functions in galaxies seems to have a 
universal powerlaw-like shape with a cut-off at luminosity of a few 
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10 40 erg s _1 and a normalization proportional to the star formation 
rate ( [Grimm et al.|2003[|L?u et al.|2006) >. This implies that most of 
the ULXs are probably a high-luminosity extension of the popula- 
tion of high-mass X-ray binaries. 

IMBHs originating from low-metallicity Population III stars 
can form binary systems due to tidal captures of single stars. The 
expected frequency of such IMBHs is not high and does not agree 
with the ULX observed frequencies (Kuranov et al. 2007 ). More- 
over, these sources should be distributed throughout the galaxy, 
while most of the ULXs are associated with the star-forming re- 
gions ( [Swartz et al.|2 009). IMBHs that are produced in a runaway 
merging in a core of young stellar clusters (Portegies Z wart et ah] 
|2004||Freitag et al.|2006] > should stay within the clusters. 

Kaaret et al. d2004| ) have found that bright X-ray sources 
(10 < Lx < 10 erg s _1 ) in three starburst galaxies are prefer- 
encially located near the star clusters. In the Antennae galaxies, 
Zeza s et aL] p002b| ) have noticed that most of the brightest X- 
ray sources are displaced from the neighbouring star clusters (they 
have considered an X-ray source to be coincident with the optical 
one if their separation is less that 2" and their astrometric accuracy 
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Figure 1. The rgb image of Antennae produced from HSTI ACS images taken in FR656N (Hen, red), F435W (blue) and F550M (green) bands. NE arrows 
length corresponds to 5". The VIMOS frames of the studied ULX regions together with Chandra positions (Of/3 radius) are indicated in green. The source 
numbers correspond to the catalogue of Zezas et al. ( 2002a). The VIMOS continuum images have been extracted in the ACS band F550M. 



was rather poor, see below). In a search for optical counterparts of 
ULXs, |Ptak et a l. (2006) have found that a few of them coincide 
with the positions of stellar clusters, but in some cases no coun- 
terparts were found. In the present paper we study the distances 
between the brightest X-ray sources and the stellar clusters in the 
Antennae galaxies which contains many ULXs. We find highly sig- 
nificant association between the ULXs and the clusters. We also es- 
timate the ages of the stellar clusters. This allows us to determine 
the minimum mass of the ULX progenitors. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-ray source sample 

The X-ray sources in the Antenna e galaxies, NGC 4038/NGC 4039 
(see Fig. Q, were surveyed by |Zezas et"aL p002a| 2006|> and 
|Swartz etaL] ( |2004| ). From the catalogue of |Zezas et al^ ( 2002a| ) we 
have selected sources situated in the main bodies of the galaxies to 
reduce the chance that they are background or foreground sources. 
We restrict the sample to those having more than 45 counts, which 
guaranties high accuracy of the coordinates. We then further down 
select the sources, which have average 0.1-10 keV luminosities in 



excess of 2.75 x 10 38 erg s 1 as estimated by|Zezas et al. (2006 1. 



These luminosities were absorption corrected assuming a power- 
law spectrum with photon index r = 1.7 and the Galactic line- 
of-sight hydrogen column density. These luminosities were recom- 
puted for the distance to the Antennae galaxies of 22 Mpc as rec- 
ommended by Schw eizer et al.|p008) . This distance gives the scale 
of 107 pc/arcsec. 

We further remove from the sample source X-37, which is 



a background quasar ( Clark et al. 2005 ) and will be further used 
for accurate astrometry. This leaves us with 20 sources (see Ta- 
ble [TJ. Sources X-24 and X-29 are situated close to the northern 
and southern nuclei, with the offsets from the bright 2MASS nu- 
clei of l'/6 and 0'/3, respectively. In addition, X-24 has a closer 
cluster counterparts. There are no indications of any AGN activity 
in the nuclei ( [Zezas et"al . 2002b, Ueda et al. 201 2]). Six sources 
with luminosities Lx > 2 x 10 39 erg s -1 ( [Zezas et al. 



2002a 



2006 1 |Swartz et aL]|2004| ) may be considered as bona fide ULXs 
(the "ULX sample"). We call other sources with luminosities in the 
range 3 x 10 38 < L x < 2 x 10 39 erg s" 1 the "sub-ULX sample". 
All the sources are located in the main bodies of the galaxies, how- 
ever, a few sources (X-18, X-19) are located in the very obscured 
dusty regions. 



2.2 Optical observations 

For spectral study we have selected regions around four ULX X-l 1, 
16, 31, 42 and also X-35. As X-32 is situated very close to X-31, 
it was also studied. We did not observe X-44, because it lies at the 
outskirts of the galaxies and we did not originally expect bright op- 
tical sources around it. Source X-29 was not observed, because it 
seemed to be associated with the nucleus of NGC 4039. We ob- 
tained two data sets on 2007 February and 2008 March with the 
Visible Multi- Object Spectrograph (VIMOS) Integral Field Unit 
(IFU) at the 8.2m Very Large Telescope Melipal at ESO's Paranal 
Observatory in Chile. 

The VIMOS IFU consists of 6400 (80x80) fibres coupled to 
microlenses and operates with four 2048 x 4096 EEV44 CCD 
detectors yielding four quadrants (Zanichelli et al. 2005 ). The 
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Figure 2. Enlarged rgb images from Fig.[T]for the regions around the brightest X-ray sources. The images were produced from ACS images taken in FR656N 
(Hen, red), F435W (blue) and F814W (green) bands. Chandra positions ($.'2 radius) of the X-ray sources are indicated in green. The bar indicates the brightest 
cluster next to an X-ray source, when it is not obvious. 
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Table 1. The brightest X-ray sources in the Antennae galaxies and the near- 
est stellar clusters. 



Object a 


RA (J2000) 


Dec (J2000) 


° a 


Offset 


Cluster d 






ULX sample 








11 


12:01:51.33 


-18:52:24.9 


39.90 


0.15 


253 










1.71 


244 










1.76 


260 


16 


12:01:52.09 


-18:51:33.2 


39.92 


1.26 


706 










1.44 


709 










1.52 


711 










1.56 


707 










1.68 


708 


29 


12:01:53.50 


-18:53:10.5 


39.57 


0.30 


49 


31 


12:01:54.28 


-18:52:01.4 


39.34 


2.50 


498 










2.61 


514 


42 


12:01:55.66 


-18:52:14.6 


39.67 


0.64 


386 










0.67 


374 










0.90 


388 










1.12 


370 










1.13 


368 










1.16 


362 


44 


12:01:56.44 


-18 51 57.4 


39.86 


3.26 


518 






sub-ULX sample 






5 


12:01:50.48 


-18:52:15.4 


38.57 


1.42 


387 










1.48 


381 










1.54 


377 


6 


12:01:50.52 


-18:52:04.2 


38.59 


0.36 


491 










0.79 


487 










1.08 


483 










1.13 


495 










1.36 


485 










1.44 


479 










1.67 


481 


13 


12:01:51.63 


-18:52:31.3 


39.09 


2.02 


206 


18 


12:01:52.40 


-18:52:06.3 


39.28 


1.92 


- 


19 


12:01:52.43 


-18:53:13.7 


38.91 


1.07 


- 


20 


12:01:52.75 


-18:51:29.5 


38.48 


0.30 


720 










0.52 


717 


24 


12:01:53.00 


-18:52:02.6 


38.80 


0.10 


500 










0.37 


501 










0.72 


503 


25 


12:01:53.01 


-18:52:09.1 


38.44 


0.67 


450 










0.88 


442 










1.40 


451 


32 


12:01:54.36 


-18:52:09.8 


38.96 


1.44 


443 










1.66 


454 










1.92 


453 


33 


12:01:54.51 


-18:53:06.3 


38.73 


0.36 


75 










0.67 


77 










1.25 


65 


34 


12:01:54.56 


-18:53:02.7 


38.90 


0.63 


90 










0.93 


89 










1.23 


87 










1.59 


88 










1.84 


86 


35 


12:01:54.78 


-18:52:51.9 


39.09 


0.57 


115 










0.83 












1.95 


109 


41 


12:01:55.49 


-18:52:35.4 


38.68 


2.24 


161 


46 


12:01:56.65 


-18:52:00.4 


38.45 


2.33 


507 



a Source number according tolZezas et al. (2002a). b Logarithm of 
luminosity in erg s — 1 in the range 0.1-10 keV jZezas et al.|2006) 
corrected for absorption and recomputed for the distance of 22 Mpc. 
c Offset (arcsec) between the X-ray source and the stellar cluster. d Cluster 
ID from the WS95 catalogue if available. 



field-of-view of \?>'! 5x\?>'! 5 was covered by 1600 spatial pixels 
(spaxels) at a spatial sampling of 0"33 per spaxel. Two different 
grisms were used, high resolution blue (^0.51 A/pixel), covering 
the wavelength range of 4150-6200 AA and high resolution orange 
(-0.60 A/pixel) covering the range of 5250-7400 AA. The spectral 
resolution measured using the sky emission lines is 1.86 and 1.88 A 
for the blue and orange grisms, respectively. We have taken spectral 
images of each target with exposure times of 6x900 s in the blue 
and 4 x 900 s in the red regions. 

For each VIMOS quadrant all the spectra were traced, identi- 
fied, bias subtracted, flat field corrected, corrected for relative fiber 
transmission, and wavelength calibrated using the routines of the 
ESO Recipe Execution pipeline (ESOrex)QWe have checked the 
fiber traces in the images and corrected them, when it was neces- 
sary, using bright night sky and nebular lines. Cosmic rays and bad 
pixels were identified and cleaned using standard IRAF[jroutines. 
Finally, the intensity of the night- sky emission line [O I] A5577 was 
used to correct for the different relative transmission of the VIMOS 
quadrants. 

The processed spectra were organized in data cubes using the 
tabulated correspondence between each fiber and its position in the 
field of view. The spectra were co- added after correcting for the 
position offset. The accuracy of the offset is ±0.2 pixel (~0 / . / 07). 
This slight deterioration of the spatial resolution does not affect the 
results. 

We use also archival imaging observations taken with the 
Hubble Space Telescope (HST). The ACS-WFC images taken in 
F435W, F550M and F814W filters on 2005 July 21 and images 
in FR656N filter taken on 2005 June 29 were used for astrome- 
try and for illustrations (Figs [T] and [2]). The WFPC2 images taken 
on 1996 January 20 in F 336W, F439W, F555W and F814W filters 
(Whit moTeTTaTl[T999t were used for photometry of the clusters 
associated with the ULXs. The WFPC2 images were reduced in 
a standard way with hstphot (Dolphin 2000). The same routine 
was used to measure all point-like sources in the images and to 
calibrate them in the standard U, B, V and / magnitudes. Using 
standard I RAF tools we carried out aperture photometry of stellar 
clusters presumably associated with the ULXs. Corrections of the 
WFPC2 magnitudes for interstellar extinction were made accord- 
ing to |Holtzman et al.| ( |1995] ). 



2.3 Astrometry of X-ray sources and stellar clusters 

To have an accurate absolute astrometry in optical images we used 
the archival r-band images taken with the MegaCam at the Canada- 
France-Hawaii Telescope (CFHT). For an astrometric solution, we 
used 15 point- like 2MASS ( |Cutri et aT]|2003| ) sources surround- 
ing the Antennae galaxies in 4' x 5' field of the CFHT images. 
The accuracy of coordinates (rms) is better than Of! 15. In the next 
step we transfered this solution from CFHT to ACS images of the 
galaxies with 22 common point-like objects, the accuracy of this 
transfer is much better than Of! 15. The same has been done with 
the WFPC2 images of the galaxies. We used a routine described 



1 ESOrex is developed and maintained by the European Southern Obser- 
vatory. 

2 I RAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy (AURA) under cooperative agreement with the National Science 
Foundation. 
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Figure 3. Positions of the bright X-ray sources (diamonds) and the stel- 
lar clusters from the WS95 catalogue (red circles). The X-ray sources are 
numbered from the catalogue o fZezas etal. ( 2002a). The size of the circle 
indicates the brightness of the cluster. The dots outline the main body of the 
galaxies. 



by[Gooch (1997). Our final accuracy in absolute astrometry of the 
optical ffiT images is the same as that in 2MASS, i.e. about 0'/2. 

The brightest 12 X-ray sources have >100 net source counts 
and the rest of the sources still have >45 cts in individual Chandra 
images. They are point-like sources and the internal accuracy of 



their coordinates is expected to be between O'/l and Of! 2 (Swartz 
|et al.|[2004] >. For the ULX- sample, we have averaged the source 
coordinates presented in surveys by Zezas et al. ( 2002a) and |Swartz| 
|et al.|pOQ4) . The ULX coordinate difference in these two surveys is 
below Of! 2 with one exception for X-35 where it is tt0f.'3. The final 
internal accuracy of the X-ray coordinates is ttOf.'l. For the sub- 
ULX-sample we adopt the coordinates from Zezas et al. (2002a). 
To determine the offset between the X-ray and optical images we 
used a background quasar ( Clark et al. 2005), the source X-37 in the 
Antennae galaxies, which is a very bright point-like source both in 
the X-ray and optical images. The coordinates of the X-ray sources 
corrected for this offset are presented in Table [T] We estimate the 
final coordinate accuracies between HSTI 'ACS and Chandra to be 
f^0 / . / 2-0 / /3, depending on the brightness of the source. 

Whitmore & Schweizer (1995 hereafter WS95) present a cat- 
alogue of stellar clusters obtained using ffiT/WFPCl with the co- 
ordinates relative to the source N442. We recomputed the coordi- 
nates of all the clusters using our solution. For all clusters in the 
vicinity of X-ray sources (see Table [TJ we have obtained indepen- 
dent measurements of the coordinates using ACS images and found 
the rms displacement from the corrected WS95 values of Of! 12. This 
allows us to use the whole WS95 catalogue with the corrected co- 
ordinates for the statistical studies and modelling of the distribution 
of the offsets between X-ray sources and clusters. 



3 ASSOCIATION WITH CLUSTERS 

In Fig. [T] we show rgb image of the Antennae galaxies with the 
brightest X-ray sources. In Fig. [2] we show enlarged rgb images 
around the X-ray sources. The positions of stellar clusters from 
WS95 catalogue together with the X-ray sources from Table [T] are 
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Figure 4. Distribution of the offsets between the positions of the X-ray 
sources and the optical positions of nearest stellar clusters (Table^. 



Table 2. Correlation between stellar clusters and X-ray sources. 
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259 


1.0 


1.8 


15 


0.09 






1.5 


4.1 


26 


0.31 






2.0 


7.2 


35 


0.63 


22.0 


371 


1.0 


2.6 


17 


0.18 






1.5 


5.8 


32 


0.5 






2.0 


10.3 


42 


1.1 


22.5 


535 


1.0 


3.7 


18 


1.2 






1.5 


8.4 


35 


0.8 






2.0 


14.9 


48 


2.0 


24.5 


690 


1.0 


4.8 


20 


2.3 






1.5 


10.8 


37 


1.2 






2.0 


19.2 


52 


2.9 



shown in Fig. [3] A remarkable property of the X-ray sources is they 
are located close to young stellar clusters and complexes. Only a 
few sources (X-ll, X-20, X-24, X-29) can be claimed to coincide 
with the position of stellar clusters (within 0"3), for other sources 
the distances to the clusters is above the astrometric accuracy. In 
Fig. [4] we show distributions of the offsets between the positions of 
the X-ray sources and the nearest stellar clusters. All the sources 
have separation below 2 / . / 7 (^290 pc), except X-44 (^350 pc). 
Source X-18 is located in a heavily obscured region and X-19 is 
at the edge of the star formation belt, with the nearest bright WS95 
clusters more that 4 /; way. However, for both these sources there 
exist much closer non-catalogued clusters (see Fig. [2}. 

The displacement between the X-ray source positions and op- 
tical clusters in Antennae was noticed by |Zezas et al.| ( [2002b| ) and 
discussed by |Zezas & Fabbiano (2002). However, they have used 
the absolute astrometry of the WFPC2 images, which has an ^ l'/5 
error, and considered X-ray sources to coincide with the optical 
sources if their separation is less than 2" . Our accurate astrometry 
does not leave any doubts that the offsets of / /5-3 // are real. 

To estimate the statistical significance of the clustering of the 
bright X-ray sources near the stellar clusters we first use the usual 
two-point correlation function, counting the number n of X-ray 
sources within an angle 5 of any cluster in a given catalogue, which 
is obtained from the WS95 catalogue. We then compare this num- 
ber to the expectation n e for the null hypothesis that the X-ray 
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Figure 5. Probabilities that the association between clusters and X-ray 
sources is random. The Poisson probabilities computed from the two-point 
correlation function for various WS95 sub-catalogues (clusters lying within 
the area identified in Fig. [3} with the cuts in visual magnitude V max =21.5, 
22.0, 22.5 and 24.5 are shown by solid, dotted, dashed and dot-dashed lines, 
respectively. 
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Figure 6. Cumulative distribution of the displacements between X-ray 
sources and stellar clusters from the WS95 catalogue (upper blue his- 
togram). The red solid histogram shows the median distribution obtained 
from 10 5 catalogues of X-ray sources distributed randomly over the galax- 
ies body within the limits shown in Fig. [3] The dotted lines show the 90% 
confidence limits. 



sources are distributed randomly over the galaxies bodies 

n e = N c i(Vm ax )N X 7r5 2 /Z, (1) 

where N c \(Vmax) is the number of clusters of magnitude below 
the cut V < Vmax, Nx is the number of X-ray sources, and 
E rsj 900CO" is the area containing stellar clusters (see Fig.[5|. We 
introduce the measure of the signal p as the probability of sampling 
n or more hits from the Poisson distribution at expectation n e . This 
probability depends on S and cuts in the stellar cluster catalogue 
VJnax- The null hypothesis is rejected with a very high significance 
p < 10 -4 for any tried 5 and any reasonable cut Vm ax (see Fig. [5] 
and Table [2]). It reaches the minimum at the 10 ~ 8 level for the cor- 
relation radius S ~0 / /7-l / /5. The deepest minimum corresponding 
to the cuts Vmax = 21.5 and 22.0 indicates that the strongest corre- 
lation exists between X-ray sources and the brightest clusters. Inter- 
estingly, that 15 out of 18 X-ray sources and all ULXs have nearest 
WS95 clusters brighter than V = 22.3. In more that half of the 
cases, the nearest WS95 cluster turned out also to be the brightest 
source with a few arcseconds. Besides that, practically all the clus- 
ters and cluster complexes, associated with the X-ray sources, are 
young (see Fig. [2]), less than 10 Myrs, because they are surrounded 
by strong Ha emission (Leitherer et al. 1999 ). This strengthens the 
association and makes the correlation revealed stronger. 

We can also check the hypothesis that the X-ray sources are 
distributed according to the stellar density (Range lov et al.|20TT] ). 
The stars, however, are distributed rather homogeneously over a 
larger area than the stellar clusters (Whitmore et al. 1999) and 
therefore the significance of the association would be even larger. 

Fig. [6] shows the cumulative distributions of the offsets be- 
tween the X-ray sources and stellar clusters from WS95 catalogue 
and those obtained from artificial X-ray source samples scattered 
randomly over the main body of the galaxies (limited by the dotted 
line in Fig.[3j. We see that the observed distribution is significantly 
higher than the median of simulated offset distributions and even 
the distribution corresponding to the 90% confidence limit. Specif- 
ically, the observed distribution has about twice as much sources 
within l'/5 than the simulated one. This confirms the association 
between bright X-ray sources and stellar clusters. 



4 CLUSTER AGES AND THE MASSES OF ULX 
PROGENITORS 

The VIMOS images (Fig.JTJ cover practically all the ULX sources 
in the Antenna galaxies. The obtained spectra were first used to 
find the extinction values from the Ha/H/3 flux ratios measured 
in the nebulae emission in the VIMOS fields around the clusters 
(Table [T] and Figs |A2f[A6| ). In the Case B of gaseous nebulae this 
ratio is 2.87, and stays the same to within 10% in a wide range 
of temperatures and densities (Osterbrock & Ferland 2006 ). Using 
Starburst99 code ( [Leitherer et al.||1999| > we calculated the evolu- 
tion tracks for the instantaneous starburst model with the Salpeter 
initial mass function for solar and double-solar metallicities. Fig. [7] 
presents the tracks together with the cluster photometric data cor- 
rected for the extinction. We then find an improved estimate for 
the extinction after matching the evolutionary tracks and the dered- 
dened V magnitudes of the clusters. For cluster N518 next to X-44 
we use the Galactic extinction Ay — 0.155 (which is the maximal 
possible value to fit the track). We note that measured Ay are in 
perfect agreement with the independent estima tes from the hydro- 
gen column density Ay — 5.5 Nu/10 22 cm 2 |Predehl & Schmitt 
1995 ) obtained from the Chandra spectra of the X-ray sources (see 
Table 5 in |Zezas et al.|2002a| >. 

We estimated cluster ages using solar metallicity as was found 
by |Bastian et al.|p009| > in their study of 16 young clusters in the 
Antennae. We find that all clusters associated with the ULXs are 
extremely young with the age of less than 5 Myr (see Table [3j. 
Using the evolutionary tracks in the (V-I)-My plane (Fig. [7]), we 
have made a horizontal shift to determine the cluster masses. 

In the X-31/X-32 VIMOS image (see Fig.|A4} there are also 
two bright stellar complexes N455 and N418 previously studied by 
Bastian et al. (2006) (named 5 and 6 there) and showing spectra 
with strong Wolf-Rayet emission features. This confirms the young 
age (< 4 Myr) of the clusters in that region. In Fig. |Al| we present 
VIMOS spectra of the clusters studied (together with spectra of 
the complexes N 5 and N 6) demonstrating very strong hydrogen, 
[O III], [Anil] and other nebular lines, which demand strong pho- 
toionizing continua of the clusters and confirm their young age. 
However, the cluster age estimates from the photometric data are 
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U-B My 

Figure 7. (a) Colour-colour and (b) colour-magnitude (for a cluster of 10 6 solar masses) diagrams of the cluster positions together with Starburst99 (Leitherer| 
|et al.|19 99i evolution tracks. The tracks are for instantaneous starburst model with the Salpeter initial mass function for solar metallicity Z = 0.02 (solid line) 
and for double-solar metallicity (dashed line). Black crosses and blue circles on the evolution tracks and numbers next to them indicate the age in million years 
(from 1 to 6 Myrs). The arrow shows a direction of the shift due to the interstellar extinction. Positions of the clusters studied are shown by small circles with 
error bars. The cluster masses shown in Table[3]were derived by horizontal shifts from the absolute magnitude My to the track position. 



Table 3. Properties of the stellar clusters next to the brightest X-ray sources in the Antennae galaxies and 
the minimum masses of ULX progenitors. 



Object" 


Cluster 6 




yd 


M v e 


M clu J 


Ages 


M h 

prog 


11 


253 


0.9-1.0 (0.93) 


18.85 


-12.86 


1.0 


4.3±0.3 


45 


16 


706 


< 0.6 (0.52) 


19.71 


-12.01 


0.5 


2.5±0.2 


100 


31 


498 


0.5±0.2 (0.62) 


18.35 


-13.36 


1.8 


2.9±0.3 


90 


32 


443 


0.8±0.2 (0.62) 


18.91 


-12.80 


1.1 


3.2±0.5 


55 


35 




3.7±0.3(3.88) 


18.33 


-13.38 


1.5 


4.6±0.4 


40 


42 


386 


0.8±0.2 (0.78) 


16.94 


-14.77 


6.7 


3.0±0.3 


75 


44 


518 


...(0.155) 


21.96 


-9.75 


0.07 


2.4±0.5 


100 



"Source number according to Zezas et al. <2002ai. Nearest stellar clusters and their ID from the WS95 catalogue if available. Extinction value determined 
from the VIMOS spectra and the corrected values (in brackets) that match the evolutionary tracks. d Dereddened WFPC2 V magnitude. e Absolute V 
magnitude for distance of 22 Mpc. -^Mass of the stellar cluster in 1O 5 M0 estimated using mass-to-light ratio M/L = 1 (in solar units). 9 Age, of the stellar 
cluster in Myr. ^Minimum mass of the ULX progenitor in Mq obtained from maximum possible age of the cluster. 



notably more accurate than those obtained from the spectroscopy. 
We note that the independent evaluation of Ay from the spectra 
break the known Ay-age degeneracy (see e.g. Bastian et al. 2009 ). 

The clear association of the X-ray sources both in the ULX 
and sub-ULX samples with the young star clusters in the Anten- 
nae galaxies indicates that these sources originate from the massive 
binaries ejected from the star clusters. The extremely young ages 
of the clusters associated with the ULXs put the lower limit on the 
mass of their progenitor stars, which varies between 40 and 100 
M (see Table[3]), indicating that all studied ULXs in the Antennae 
galaxies are associated with the most massive stars. 



5 DISCUSSION 

The correlation between X-ray source positions and stellar clusters 
was already noticed by Kaaret et al. ( 2004 ) for three starburst galax- 
ies M 82, NGC 1569 and NGC 5253 (which lie closer than 4 Mpc). 
In spite of their astrometric accuracy of >1 ;/ , the proximity of the 



galaxies permitted them also to find the displacement between X- 
ray sources (with X-ray luminosities less than 10 39 erg s _1 ) and 
bright clusters. They found that brighter X-ray sources preferen- 
tially occur closer to the clusters, and no sources at luminosities 
above 10 38 erg s _1 are displaced by more than 200 pc. The lower 
luminosity sources are mostly low-mass X-ray binaries, which to- 
tally loose connection to their parent clusters because of their large 
ages or because the clusters disappear on that time frame. In their 
study only one galaxy (M 82) has sources brighter than 10 38 erg 
s _1 and only one source was brighter than 10 39 erg s _1 . As the 
Antennae galaxies have a relatively high star-formation rate, most 
of bright X-ray sources studied by us are probably high-mass X-ray 
binaries. 

Recently Rang elov et al.| ( |2011|> found that bright sources 
(with luminosities 3 x 10^-8 x 10 3 * erg s" 1 ) in NGC 4449 are 
located close to the young stellar clusters with the age less than 6-8 
Myr. This association was, however, only marginal, when they as- 
sumed that X-ray sources follow the light of the galaxy. It is worth 
mentioning that about half of the X-ray sources in NGC 4449 are 
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probably low-mass X-ray binaries and the correlation between the 
remaining candidate high-mass X-ray binaries and stellar clusters 
is significant. 

In the Antennae galaxies, Clark et al. (201 1 ) identified 38 in- 
frared counterparts within the l'.'O error boxes of 120 X-ray sources 
from the Zeza s et ah] {2006) catalogue. The majority of the IR 
sources are star clusters. For ten of them they were able to measure 
masses of about 1O 6 M and ages, which lie in a very narrow range 
around 10 7 yr. However, their study covers only clusters within the 
error circles, which dramatically reduces the statistics for the bright 
X-ray sources. 

The recent study by Rangelov et al. (2012) showed that 22 out 
of 82 X-ray sources from the same catalogue have star clusters as 
optical counterparts within 1" 1 error box[jThey have measured the 
age of these clusters and the extinction using UBVI and Ha mag- 
nitudes finding that 14 sources are situated next to young clusters 
with ages of 4-6 Myr. Among the clusters associated with ULXs, 
they found the age of 5 Myr for X-ll (their source X27), which 
is in agreement with our estimate of 4.3 =b 0.3 Myr. However, for 
X-42 (their X99), they find the age of 200 Myr, while we find 3 
Myr. This discrepancy probably results from the Ay-age degen- 
eracy in colour-colour diagrams (see e.g. |Bastian et al.|2006| ). On 
the other hand, our measurements of the Ay are independent of the 
photometry. 

The study of Rang elov et al.| < |2Q12) > is limited only to the 
clusters very close to the X-ray sources. We show, however, that 
even clusters that are further away are also associated with the X- 
ray sources and their conclusion that black holes stay within the 
clusters contradicts our results. We find that majority of the X-ray 
sources with luminosities > 3 x 10 38 erg s _1 are located close 
to (but not within) very young clusters. Together with findings by 
Kaar et et"ah] ( |2004| ) this means that the association with clusters 
and the displacements are the common properties of all high-mass 
X-ray binaries and the ULXs do not differ in this aspect from the 
less bright X-ray binaries. The correlation and the displacement of 
the brightest X-ray sources up to 300 pc from the clusters imply 
that they have been ejected from the clusters. It is not possible that 
the X-ray sources have been left behind their dissolved host clus- 
ters, because the extremely small ages of the clusters do not leave 
any time for the cluster (and nearby clusters) disintegration. There 
exist three known runaway scenarios, which can explain these dis- 
placements: 

(i) Instantaneous, symmetric mass loss in a supernova explosion in 
a close binary (|Zwicky|1957} |Blaauw|1961| [Shklovskii |1976| |van| 
|den Heuvel et d.|200fl^ " 

(ii) Instantaneous momentum impulse, or "kick", imparted through 
asymmetry of the supernova explosion in formation of a neutron 
star ( |Woosley|1987l|Cordes & Chernoff|1998] ). A non-radial asym- 
metry in neutrino momentum distribution of only 1 % would give 
the neutron star a recoil velocity of a few hundred kilometers per 
second. 

(iii) The ejection of massive X-ray binaries from the clusters due 
to close 3- and 4-b ody encounters ([Poveda et al.||1967| [R eggie 
T975l|Mikkola|1983HMoeckel & Bate|2010||Pflamm-Altenbu rg & 



Kroupa|2010| L 



3 Their astrometry has la error of Of/ 54, but we note that the difference 
between the HST and Chandra coordinates of the QSO X-37 is certainly 
much smaller than quoted value of the mean shift of 2'/ 5 175 obtained by 
[Rangelov et aL|(2012} . 



The first two mechanisms cannot explain the data, because the 
binary is ejected after the SN explosion, and the system stays in 
the cluster for all the period of the primary evolution. This time is 
« 3 Myr for « 100M© star ( [Leitherer et al.|1999[ >, which leaves 
no time for relocation to their present positions. The direct SN kick 
(ii) deals with neutron stars, but even rescaling the kick velocity 
relation to a black hole of ~ 1OM ( Zezas & Fabbiano 2002) gives 
velocity of less than 3 km s _1 , because the companion, which also 
had to evolve, has to be very massive. 

Alternatively, the ejection of massive X-ray binaries from the 
clusters can occur already at the very beginning stages of the clus- 
ter evolution due to close encounters. This scenario is supported 
by the presence of the blue halos of massive stars around the clus- 
ters ( Whit more et al. 1999 ) as seen in Fig[2] The belts of young blue 
stars close to the young clusters might indicate that some of the less 
massive clusters have already dissolved (Moeckel & Bate 2010| ). 
Due to gas expulsion, young embedded star clusters become super- 
virial and start to expand on the typical time- scale of « 1 Myr 
( |Kroupa et al.|2001|[Prlamm-Altenburg & Kroupa|2 006 ). The rela- 
tively large displacement and high spatial velocities are consistent 
with the results of N-body simulations showing that the most mas- 
sive stars (binaries) are ejected at early stages of formation of stellar 
clusters ( [Moeckel & Bate|2010| ). 

Among the ULX sample (Table [3} only one source X-ll re- 
sides in a cluster. Other six sources have an average ejection ve- 
locity (accounting for projection) of 83 km s _1 . Five of these six 
sources are associated with the extremely young clusters of ages 
from 2.4 to 3.2 Myr, with the mean age of 2.8 Myr (corresponding 
to the progenitor masses of > 100M Q ). The average displacement 
from the clusters is l'/94, giving the deprojected mean distance of 
260 pc and the mean ejection velocity of 92 km s _1 . 

The short cluster life-times imply limitation on the ULX bi- 
nary mass ratio. The main sequence time for a ~ 1OOM star is 
w 3 Myr. One may understand the association of the five ULXs 
with young clusters (see Table[3]), if both companions are very mas- 
sive Mi r>j M 2 ~ 1OOM and they evolve at approximately the 
same time-scale T\. If the original mass ratio is large Mi ^> M2, 
the evolution of the pre-ULX system takes twice as long ~ 2Ti , 
because after time T± the primary transfers its mass to the com- 
panion. Depending on the original mass ratio (and the masses) the 
total binary evolution takes between T\ and ~ 2Ti . Two sources, 
X-ll and X-35, may be associated with less massive binaries with 
Mi -5O-1OOM and M 2 - 5OM . 

The observations, however, are clearly inconsistent with the 
IMBH scenario as such massive black hole cannot be ejected from 
clusters at the required high (80-90 km s _1 ) velocities. In their dis- 
placements from young clusters, the ULXs do not differ from the 
sub-ULXs sources. In this property, the ULXs constitute a contin- 
uation of the high-mass X-ray binary population. 



6 CONCLUSIONS 

In this paper we have studied the environments around the bright- 
est X-ray sources in the Antennae galaxies. Using a high-accuracy 
astrometric solution we found a highly significant association be- 
tween the X-ray sources and the stellar clusters. We also showed 
that most of the bright X-ray sources are located outside of these 
clusters confirming previous findings. We have studied two samples 
of the X-ray sources, the ULX-sample (L x > 2 x 10 39 erg s -1 ) 
and sub-ULX one (3 x 10 38 < L x < 2 x 10 39 erg s" 1 ). Using 
VIMOS imaging we were able to determine the reddening towards 
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to the clusters near the ULXs and to find their age, which was less 
than 5 Myr. The sub-ULX sample sources are also connected with 
young clusters (< 10 Myr) because they surrounded by Ha emis- 
sion. 

The clear association of the young star clusters and the X-ray 
sources both in the ULX and sub-ULX samples in the Antennae de- 
mands these massive binaries to be ejected off the star clusters. All 
these sources belong to the runaway star population. The extremely 
young ages of the clusters associated with the ULXs put the lower 
limit on the mass of their progenitor stars. The minimum possi- 
ble mass varies between 40 and 1OOM , indicating that all studied 
ULXs in the Antennae galaxies are associated with the most mas- 
sive binary stars. A few ULX sources associated with the youngest 
star clusters must originate from ~ 1OOM stars with the mass ra- 
tio ~ 1. We conclude that the association with young star clusters 
and the displacements are the common properties of all high-mass 
X-ray binaries and the ULXs do not differ in this aspect from the 
less bright X-ray binaries. Our findings strongly rule out an alter- 
native hypothesis that majority of ULXs are IMBHs. 

The star ejection mechanisms based on a SN explosion (either 
symmetric explosion in a close binary or a direct "kick", imparted 
through asymmetry of the SN) are ruled out because the binary is 
ejected after the SN explosion. Because the ULX-clusters are very 
young, 3-4 Myr, there is no time for the pre-SN evolution inside 
the cluster. 

The relatively large displacements (up to 300 pc) and high spa- 
tial velocities (80-90 km s _1 ) found are consistent with the results 
of N-body simulations showing that the most massive stars (bina- 
ries) are ejected at the very beginning stages of the cluster evolution 
due to close encounters. 
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Figure Al. The VIMOS spectra of the clusters associated with X-ll, 16, 31, 32, 35 and 42 (see Table[3j as well as the cluster complexes N455 and N418 
(objects 5 and 6 in Bastian et al. 2006 see Fig. |A4) taken with the blue (panel a) and orange (panel b) grisms. The main spectral lines are indicated. 



APPENDIX A: VIMOS IMAGING, SPECTRA OF 
CLUSTERS AND COMMENTS ON INDIVIDUAL 
SOURCES 

Fig. |A1| presents the spectra of the bright stellar clusters next to 
the selected X-ray sources. In Figs |A2f[A6] we show VIMOS im- 
ages of regions around six brightest X-ray sources. The images 
were derived in Ha, [S II] A6731, [Om] A5007. There are also 
extinction Ay maps and images in the continuum extracted in the 
ACS/F550M band and the ACS F550M image itself for compar- 



ison. In the estimates of the reddening we used the Ka/H/3 flux 
ratios measured from the nebulae emission in the VIMOS fields. 
We adopt Rv = 3.1 in the estimates. 

The source X-18 (see Fig. [2]) is located in a dusty zone, where 
there are no clusters from WS95 catalogue. There are also no 
known clusters around the source X-19. Source X-35 is located in 
a dust lane with two clusters at about the same distance (0 / /6-0'. / 8) 
nearby. We have studied (Table [3} the brighter object, which is not 
from the WS95 catalogue, because another cluster is too weak for 
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Figure A2. VIMOS images of the region around X-ll with the X-ray position marked by circles (with 0'/5 radius). The images show the line intensities in 
Ha, [S II] A6731, [O in] A5007, the extinction A v , the continuum flux measured by VIMOS in the ACS F550M band and the ACS F550M continuum image 
itself. The line intensities (measured in 10 _ 16 erg cm -2 s _1 pixel -1 ) were not dereddened. North is to the right and east is up. 
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Figure A3. Same as Fig.|A2| but for the source X-16. 
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Figure A4. Same as Fig. |A2| but for the sources X-31 and X-32. Two stellar complexes N455 and N418 (objects 5 and 6 from the study of Bastian et al.|2006) 
are located ~3 x /5 to the NE and ~5 /; to the SE from the X-ray positions of X-32, respectively. The nearest clusters to X-31 and X-32 are situated ~3"SW and 
~1 /; NE from the corresponding source. 
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Figure A6. Same as Fig.|A2| but for the source X-42. 



a detailed study. The objects associated with the X-ray sources X- 
16 and X-42 are cluster complexes and the brightest cluster in the 
complex was selected. 



Al X-ll 

This X-ray source is located inside a cluster with dereddened 
V = 18.85 (see Figs |2| and |A2[ Table [TJ. The cluster's spectrum 
(see Fig. \A\\ has strong emission lines, which are formed in H II 
regions, however, Ha line has broad wings indicating their stellar 
origin. In addition to the hydrogen lines Ha and H/3 we see lines 
[Ol] AA6300, 6363, [Om] AA4959, 5007, [Sn] AA6716, 6730, 
[Nil] AA6548, 6583, [Nil] A5754. There are also lines indicating 
strong UV ionization sources, they are He I AA4922, 5015, 5876, 
6668, 7065, [Ni] AA5197, 5200 and [Ann] A7135. Although the 
[O I] lines are not resolved, other nebular lines are slightly broad- 
ened. Correcting for the spectral resolution we find FWHMs of 
hydrogen, [Om], [Arm] lines of 50, 70 and 80km s -1 respec- 
tively. There are numerous Fe II narrow lines both in emission and 
in absorption and absorption line of Mgll A4481. The strongest 
Fell A5169 line shows clear P Cygni-like profile. The absorption 
lines indicate their origin in stellar atmospheres. In the red spectral 
region, the Fe II lines appear in emission. We suspect a presence 
of the [Fell], Nil, Nm and Sll emission lines in the spectrum, 
however, the spectrum quality is not enough to study these lines. 



A2 X-16 

This X-ray source is located close (l'/5) to a faint cluster with 
dereddened V = 19.71 (see Fig s |2| and | A3 [ TablefTJ. The nebular 
lines in the cluster spectrum are about the same as those in cluster 
next to X-l 1, but the [O III] A5007 line is brighter that H/3, which 
indicates a stronger UV ionization source. The narrow Ha shows 
broad wings. The Fell absorption lines and weak Sill AA6347, 
6371 lines are present. 



A3 X-31 and X32 

The spectra of the clusters located close to the X-ray sources X-31 
(N498, dereddened V = 18.35, separation 2 / /5) and X-32 (N443, 
V = 18.91, separation l'/5) are shown in Fig. Al together with 
th e spectra of two othe r clusters N455 and N418 (named 5 and 6 
in |Bastian et al.|2006| ). Cluster N455 is located - 3" 5 to the NE 
and the brigth cluster N418 is located ~ 5"to SE from X-ray posi- 



tions of X-32 (see Figs[l)|2|and |A4| l. We show spectra of these four 
clusters together because they are similar and all clusters are lo- 
cated close to each other. N455 and N418 show bright Wolf-Rayet 
spectral features, He II A4686, Bowen blend (C m/N III emission 
lines, where Nm lines dominate strongly) together with the nu- 
merous [Fe III] and [Fe II] emission lines. Emission lines of Fe II 
and Si III are also detected. There are N Il-N IV emission lines as 
well, however, the Wolf-Rayet hump at ~ 5800A (C IV and N IV 
lines) is not strong. The Ha lines show broad emission wings in the 
spectra. All these features point towards WN and LBV-like stars 
with extended winds in N455 and N418 clusters, indicating their 
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ages to be less than 6 Myrs. The spectra of the clusters, which we 
accosiate with X-31 and X-32, show similar features. The cluster 
next to X-31 shows Ha line with broad wings, the N III lines in the 
Bowen blend, many [Fe III], [Fe II] and N II emission lines, as well 
as [Ar IV] AA7170, 7237, 7262, 7331 lines in addition to the bright 
[Arm] A7135 line. Broad He II A4686 is marginally detected. The 
cluster next to X-32 shows similar features, however, its [Fein] 
emission lines are notably broadened, and no [Ariv] lines were 
detected in the spectrum. Emission lines of Fe II, N II and Si III to- 
gether with a faint He II A4686 hump are detected. 



A4 X-35 

This X-ray source is located close (07 6-078) to two clusters. The 
closest cluster N115 from WS95 catalogue (0 / /57 to the north) is 
too weak for a detailed study. Another apparently rather faint clus- 
ter (not m entio ned in WS95, with dereddened V = 18.33; see 
Figs [2] and [A5] Table [1} is situated 0783 to the SSW from X-35 



at the edge of a strong dust line. The cluster is notably reddened 
and its spectrum is noisy. We find about the same nebular spectrum 
as that in X-16 with the [O III] A5007 line being brighter that H/3. 
However, in the red spectral region He I nebular lines and [Ar III] 
are strong and broadened. Interestingly, there are [Fe III] emission 
lines in the spectrum like those we observe in N 5 and N 6 clusters 
close to X-32. The [Fe II] emission lines are marginally detected. 



A5 X-42 

This X-ray source in located close (at 078) to a bright cluster com- 
plex with total (dereddened) V = 16.94 (see Figs [2] and | A6[ Ta- 
ble^. The set of nebular lines is the same as in other clusters. How- 
ever, this complex includes stars of different ages, there is a notable 
amount of A-type stars, because the H/3 line has broad absorption 
wings. There are many stellar absorption lines in the spectrum indi- 
cating the B/A populations, numerous He I, Fe II, C II lines as well 
as the Mg II and Si II lines. He I lines appear in emission in the red 
part of the spectrum. 



